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Within the t-J model, the heat transport of the underdoped cuprates is studied based on the 
fermion-spin theory. It is shown that at low temperatures the energy dependence of the thermal 
conductivity spectrum consists of two bands. The higher-energy band shows a weak peak, while the 
low-energy peak is located at a finite energy. This high-energy broad band is severely suppressed with 
increasing temperatures, and vanishes at higher temperature. It is also shown that the temperature 
dependence of the thermal conductivity increases monotonously with increasing temperatures, in 
agreement with experiments. 
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It has become clear in the past fifteen years that the 
cuprate superconductors are indeed fundamentally differ- 
ent from other conventional metals in that they are doped 
Mott insulators^'^. The undoped cuprates are Mott insu- 
lators with the antiferromagnetic (AF) long-range-order 
(AFLRO). A small amount of carrier dopings to this Mott 
insulating state drives the metal-insulator transition and 
directly results in the superconducting transition at low 
temperatures for low carrier dopings^'^. The optical 
conductivity in the normal state above superconducting 
transition temperature shows a non-Drude behavior at 
low energies, and is carried by x holes, with x is the hole 
doping concentration, while the resistivity exhibits a lin- 
ear temperature behavior over a wide range of temper- 
atures in the underdoped regime^. This unusual charge 
transport does not fit in the conventional Fermi-liquid 
theory'*. Furthermore, a clear departure from the uni- 
versal Wiedemann-Franz law for the typical Fermi-liquid 
behavior is observed in doped cuprates^. It has been ar- 
gued that these anomalous features may be interpreted 
within the framework of the charge-spin separation"*'^'^, 
where the electron is separated into a neutral spinon and 
charged holon, therefore the basic excitations of doped 
cuprates are not fermionic quasiparticles as in other con- 
ventional metals with charge, spin and heat all carried 
by one and the same particles. 

The heat transport, as manifested by the thermal con- 
ductivity, is one of the basic transport properties that 
provides a wealth of useful informations on the carriers 
and phonons as well as their scattering processes*~^°. In 
the conventional metals, the thermal conductivity con- 
tains both contributions from carriers and phonons®. The 
phonon contribution to the thermal conductivity is al- 
ways present in the conventional metals, while the mag- 
nitude of the carrier contribution depends on the type of 
material because it is directly proportional to the free 
carrier density. In particular, the free carrier density 
in conventional superconducting materials is very high, 
then the thermal conductivity from the free carriers is 
usually the dominant contribution. In the cuprate mate- 
rials, the physical properties in the undoped case are now 
quite well understood* : here the system of interacting lo- 



calized Cu^+ spins is well described as the conventional 
insulating antiferromagnet, and at low temperatures, it 
has the purely phonon thermal conductivity^. However, 
it has been shown from experiments^'** that the phonon 
contribution to the thermal conductivity is strongly sup- 
pressed in the underdoped regime, and the conventional 
models of phonon heat transport based on phonon-defect 
scattering or conventional phonon-clcctron scattering fail 
to explain the experimental data**. Since the electron 
has been separated as the holon and spinon, and therefore 
it has been argued that the contributions from spinous 
and holons may dominate the heat transport of the un- 
derdoped cuprates**'''. Recently, we*^ have developed a 
fermion-spin theory based on the charge-spin separation 
to study the physical properties of doped cuprates, where 
the electron operator is decoupled as the gauge invariant 
dressed holon and spinon. Within this theory, we have 
discussed charge transport and spin response of the un- 
derdoped cuprates and superconducting mechanism. It 
has been shown that the charge transport is mainly gov- 
erned by the scattering from the dressed holons due to 
the dressed spinon fluctuation, while the scattering from 
the dressed spinons due to the dressed holon fluctuation 
dominates the spin response*^'*'^. These dressed holons 
interact occurring directly through the kinetic energy by 
exchanging the dressed spinon excitations, leading to a 
net attractive force between the dressed holons, then 
the electron Cooper pairs originating from the dressed 
holon pairing state are due to the charge-spin recombina- 
tion, and their condensation reveals the superconducting 
ground-state, where the electron superconducting tran- 
sition temperature is determined by the dressed holon 
pair transition temperature, and is proportional to the 
hole doping concentration in the underdoped regime*'*, 
in agreement with the experiments. In this paper, we 
apply this successful approach to discuss the heat trans- 
port of the underdoped cuprates. Within the t-J model, 
we show that although both dressed holons and spinons 
are responsible for the heat transport of the underdoped 
cuprates, the contribution from the dressed spinons dom- 
inates the thermal conductivity. 

In the doped cuprates, the single common feature is 
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the presence of the two-dimensional (2D) Cu02 plane-'^, 
then it is behoved that the unusual physical properties 
are closely related to the doped Cu02 planes. It has been 
argued^ that the essential physics of the doped Cu02 
planes is contained in the 2D t-J model, 

^ = E ^l^i+v- + E 

+ J Sj • Si+^, (1) 
irj 

with fj = ±i, zby, Cj^ (Cia) is the electron creation (an- 
nihilation) operator, Sj = CjaCi/2 is spin operator with 
a = {(Jx,<Jy,<Jz) as Pauli matrices, and is the chemi- 
cal potential. The t-J model (1) is supplemented by the 
single occupancy local constraint Y^^Cl^Ci^ < 1. This 
local constraint reflects the strong electron correlation in 
the doped Mott insulator^, and can be treated properly 
in analytical form within the fermion-spin theory^^ based 
on the charge-spin separation, 
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where the spinful fermion operator hia = 

describes the charge degree of freedom together with 
some effects of the spinon configuration rearrangements 
due to the presence of the hole itself (dressed holon), 
while the spin operator Si describes the spin degree of 
freedom (dressed spinon), then the electron on-site lo- 
cal constraint for the single occupancy, J2(T^ia^i<^ — 
S+h,^hl^S- + S-h^ih^S^ = heists- + SrS+) = 
is satisfied in analytical calculations, 
and the double spinful fermion occupancy, h\^h\_^ = 

e**- /it /lie***- = 0, hiahi-^ = e-'^^/ii/iie"'*'- = 0, 
arc ruled out automatically. It has been shown that these 
dressed holon and spinon are gauge invariant, and in this 
sense, they are real and can be interpreted as the physi- 
cal excitations"'^^'-'^^. This dressed holon hia- is a spinless 
fermion hi incorporated a spinon cloud e^***" (magnetic 
flux), then is a magnetic dressing. In other words, the 
gauge invariant dressed holon carries some spinon mes- 
sages, i.e., it shares its nontrivial spinon environment due 
to the presence of the holon itself^®. Although in com- 
mon sense hia- is not a real spinful fermion, it behaves 
like a spinful fermion. In this fermion-spin representa- 
tion, the low-energy behavior of the t-J model (1) can be 
expressed as^^, 

H = Y,Hi, (3a) 

i 

Hi = -tY^ih^^Sfhl^^S-^^ + /iaSr/it^.^5+ .) 

v 

- ^,^h\^hia + JeS^^i ■ ^i+fi, (3b) 

with Jeff = (1 — xYJ, and x = {h\^hia) = {h\hi) is the 
hole doping concentration. As a consequence, the kinetic 



energy (i) term in the t-J model has been expressed as 
the dressed holon-spinon interaction, which dominates 
the essential physics of doped cuprates, while the mag- 
netic energy (J) term is only to form an adequate spinon 
configuration. 

Within the t-J model (3), the thermal conductivity of 
doped cuprates can be expressed as^'', 

1 lmRQ{uj,T) 



k{(j,T) 



(4) 
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with I1q{uj,T) is the heat current-current correlation 
function, and is defined as, 
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-{TrjQ{r)jQ{r')), 



(5) 



where r and r' are the imaginary times, Tr is the r or- 
der operator, while the heat current density is obtained 
within the Hamiltonian (3b) by using Heisenberg's equa- 
tion of motion as^^, 
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(6a) 



Jq^ = i{xtf 5Z 'nh\+f^,ahi+fia + itJ-Xt^fih\^-^hi„, (6b) 
ij) = ii(eJeff)2 ^(,7 - n')[StSi_^^^,St_^ 

if}fi' 

ifjf}' 

-5+Sr.5t^+^.], (6c) 

where is lattice site, e — 1 -\- 2t<j)/Jcs, x = {S^S~^f^) 

is the spinon correlation function, and ^ = {hl^hi+fj^) 
is the dressed holon's particle- hole parameter. Although 
the heat current density jg has been separated into two 

parts jg'^ and with jg'-' is the holon heat current 

(s) 

density, and jg is the spinon heat current density, the 
strong correlation between dressed holons and spinous 
still is considered self-consistently through the dressed 
spinon's order parameters entering in the dressed holon's 
propagator, and the dressed holon's order parameters en- 
tering in the dressed spinon's propagator. In this case, 
the heat current-current correlation function (5) can be 
calculated in terms of the full dressed holon and spinon 
Green's functions ga-(k,uj) and D{k,uj). Following the 
discussions of the charge transport^^'^"^'^^, we obtain the 
thermal conductivity of doped cuprates as, 

k{lo,T) = Kh{c^,T) + Ks{uj,T), (7a) 

} J —oo 
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.Mk,.' f^^^' + ^^l-''-^^'\ (7c) 
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where K,h{LO,T) and Ks{<-^,T) are the corresponding con- 
tributions from dressed holons and spinons, respectively, 
7^^^ = (sin^fc^ + sin2fcj;)/4, A;, = Zxtifi - Zxtjk), As = 
(ZJeff)2e(2ex+2C-4x7fe), 7k = {l/Z) Y.f, e*''"*. ^ is the 
number of the nearest neighbor sites, the spinon correla- 
tion function C = {l/Z'^)Y,~ ~,{Sf^~S-^^}), the dressed 
holon and spinon spectral functions are obtained as 
Aha(k,ijj) = —21mga{k,tjj) and As{k,uj) = — 2Im£'(fc, w), 
respectively, while the full dressed holon and spinon 
Green's functions have been discussed in detail in Ref.-'^^, 
and can be expressed as g^^{k,uj) = 5^"' ^{k^uj) — 
T.h(k,uj) and D-^{k,uj) = D^°^-^{k,u;) - T.s{k,uj), with 
the mean-field (MF) dressed holon and spinon Green's 



functions. 
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- Cfc and D^°^-\k,Lo) 



(w — u)l)/Bk, and the second-order dressed holon and 
spinon self-energies are obtained by the loop expansion 

to the second-order^^'^^ as. 



p+p' 



+ 



(jj + UJp+p' — LOp' — ^p+k W + LOpi — LOp+pi — £,p+k 

Ft\k,p,p') Fi''\k,p,p') \ 



U) + LOpi + LOp^pi — ^p+k ~ '^p+p' ~ <^p' ~ ^p+k I 

(8a) 

^ ^ pp' 



Ft\k,p,p') 



F^'\k,p,p') 



\ljJ + ^p+p' — ^p' — IJ^k+p <^ + ^p+p' ~ ip' + '^k+p I 

(8b) 

respectively, where Bk = A[2x^(e7k - 1) + x(7k - e)], 
A = 2ZJeff, the spinon correlation function — 

{SfSt^-), F['\k,p,p') = nf(W)[l - riFiip')] - 
nB{u)k+p)[nF{^p') 

npiCp+p')], F!f\k,p,p') = nF{ip+p')[l - ripi^p')] + 
[1 + nB{uJk+p)\[nF{ip') - nF{S,p+p')], F['^\k,p,p') = 
Uf {ip+k)[nB (wp' ) - ns (i^p+p' )]+nB [ojp+p' ) [1 + r^s (wp' )] , 
F2^\k,p,p') = nFi^p+k)[nB{oJp'+p) - nB{(^p')] + 
nB{oJp')[l + nB{u)p'+p)], F^''\k,p,p') = nF{^p+k)[i + 
nB{ujp+p')+nB{oJp')]+nBiujp,)nB{ojp+p'), F^'^\k,p,p') = 
"F(^p+fc)[l + nB(uJp+p') + nB(wp')] - [1 + n_B(wp')][l -|- 
nB{i^p+p')], TiBi^Jp) and nF{£,p) are the boson and 
fermion distribution functions, respectively, and the MF 
dressed holon and spinon spectra are given by ^k = 
Ztxjk - Mi and luI = ^i(7fe)^ + ^27*: + A^, respectively, 
with Ai = aeA2(ex^ + x/2), ^2 = -eA2[a(x" + ex/2) + 
{aC- + (1 - a)/{AZ) - a€x/{2Z)) + {aC+{l - a)/{2Z) - 
ax' 12)12], A3 = \^[aC' + (1 - a)/(4Z) - aex/(2Z) + 
e^{aC a)/{2Z) - ax72)/2], and the spinon corre- 

lation function = (1/^^) E^,v(^^>5^t^-')- o^^^er 



not to violate the sum rule of the correlation function 
{S^S') = 1/2 in the case without AFLRO, the impor- 
tant decoupling parameter a has been introduced in the 
MF calculation^^'^°, which can be regarded as the vertex 
correction. All the above MF order parameters, decou- 
pling parameter q, and chemical potential /U are deter- 
mined by the self-consistent calculation^". 

The thermal conductivity is one of the direct probes 
to observe the low energy quasi-particles through its fre- 
quency and temperature dependences^"^". In Fig. 1, we 
present the results of the thermal conductivity k{u)) as 
a function of frequency at doping x = 0.06 (solid line), 
X = 0.10 (dashed line), and x = 0.12 (dash-dotted line) 
for parameter t/J = 2.5 with temperature T = 0.05J. 
Although k{uj) is not observable from experiments, its 
features will have observable implications on the observ- 
able k{T). From Fig. 1, we find the thermal conductiv- 
ity spectrum consists of two bands separated at a; ~ 0.5t, 
the higher-energy band, corresponding to the midinfrared 
band in the optical conductivity, shows a weak peak at 
a; ~ It = 2.5J, while the position of the lower-energy 
peak in the present underdoped cuprates is doping de- 
pendent, and is located at a finite energy lu ^ xJ. More- 
over, we also find from the above calculations that al- 
though both dressed holons and spinons are responsible 
for the thermal conductivity k{lv), the contribution from 
the dressed spinons is much larger than these from the 
dressed holons, i.e., Ks{uj) ^ Kh{oj) in the underdoped 
regime, and therefore the thermal conductivity of the 
underdoped cuprates is mainly determined by its dressed 
spinon part Ks{lu). For a better understanding of the 
heat transport of the underdoped cuprates, we have stud- 
ied the frequency dependence of the thermal conductiv- 
ity spectrum at different temperatures, and the results 
at a; = 0.10 for t/J = 2.5 in T = 0.05J (solid line), 
T = 0.25 J (dashed line), and T = 0.5 J (dash-dotted 
line) are plotted in Fig. 2. These results show that the 
high-energy band in k{ui) is severely suppressed with in- 
creasing temperatures, and vanishes at higher tempera- 
ture (T > 0.4 J). 

Now we turn to discuss the temperature dependence 
of the thermal conductivity k{T), which can be obtained 
from Eq. (7) as k{T) = lim^^o k{lo,T). The results 
of k{T) a,t X = 0.10 (solid line), x = 0.12 (dashed 
line), and x = 0.15 (dash-dotted line) for t/J = 2.5 
are shown in Fig. 3 in comparison with the experimen- 
tal results'^ taken on La2-xSr2:Cu04 (inset). Our re- 
sults show that the thermal conductivity k{T) increases 
monotonously with increasing temperatures for T < xJ, 
and is very weak temperature dependent for T > xJ, 
in good agreement with the experimental data^ in the 
normal state, where there is a shoulder around the tem- 
perature T ^ xJ for both experimental and theoretical 
results, and is consistent with the position of the lower- 
energy peak in k{u)). On the other hand, we have noted 
that the smooth evolution of k(T) from the supercon- 
ducting state (T < 0.02J « 20K) to the normal state 
(T > 0.02 J) has been observed from the experiment^ for 
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La2-a;Sr2;Cu04- Based on the superconducting mecha- 
nism driven by the kinetic energy the thermal conduc- 
tivity in the superconducting state {T < 0.02 J w 20K) 
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FIG. 1. The thermal conductivity as a function of fre- 
quency ai X = 0.06 (solid line), x = 0.10 (dashed line), and 
X = 0.12 (dotted line) with t/J = 2.5 in T = 0.05J. 
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FIG. 2. The thermal conductivity as a function of fre- 
quency at a; = 0.10 in r = 0.05 J (solid line), T = 0.25 J 
(dashed line), and T = 0.5 J (dotted line) with t/J = 2.5. 
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FIG. 3. The thermal conductivity as a function of temper- 
ature a,t X = 0.10 (solid line), x = 0.12 (dashed line), and 
X = 0.15 (dotted line) with t/J = 2.5. Inset: the experimen- 
tal result of La2-2,Sr2,Cu04 taken from Ref. [9]. 



is under investigation now. 

In the above discussions, the central concern of the 
thermal conductivity in the underdoped cuprates is the 
charge-spin separation, then the heat transport is mainly 
determined by the contribution from dressed spinous 
Ks{uj,T). Since K,a{^iT) in Eq. (7c) is obtained in 
terms of the dressed spinon Green's function D{k,uj), 
while this dressed spinon Green's function is evaluated 
by considering the second-order correction due to the 
dressed holon pair bubble^^'^^, therefore the observed un- 
usual frequency and temperature dependence of the ther- 
mal conductivity spectrum of the underdoped cuprates is 
closely related to the incommensurate spin dynamics^^. 
This is because that within the fermion-spin theory, the 
dynamical spin structure factor has been obtained^^'^^ in 
terms of the full dressed spinon Green's function as, 



S'(k,cj) = -2[l + nBiuj)]lmDik,uj) 
= [l + nBiL^)]Asik,uj) 

^ -2[1 + nB{uj)]BkImJ:s{K^) 

[w2 ^u;l- Rcl],(k,c^)]2 + [ImE,(k,c^)] = 



(9) 



where ImEs(k,a;) and ReSs(k, w) are corresponding 
imaginary part and real part of the dressed spinon self- 
energy function Ss(k, w) in Eq. (8b). As we have 
shown in detail in Refs.^^'^^, the dynamical spin struc- 
ture factor (9) has a well-defined resonance character. 
5(k, a;) exhibits a peak when the incoming neutron en- 
ergy uj is equal to the renormalized spin excitation E"^ — 
ujI + BkKel^s{k, Ek) for certain critical wave vectors k^ 
(positions of the incommensurate peaks). The height 
of these peaks is determined by the imaginary part of 
the dressed spinon self-energy 1/Iml]s(k5, w), i.e., the 
height of the incommensurate peaks is determined by 
damping, it then is fully understandable that they are 
suppressed as the energy and temperature are increased. 
Since this incoming neutron resonance energy u = Ek is 
finite, this leads to that the lower-energy peak in k(w) 
is located at a finite energy. Near the half-filling, the 
spin excitations are centered around the AF wave vec- 
tor [7r,7r], so the commensurate AF peak appears there. 
Upon doping, the dressed holons disturb the AF back- 
ground. Within the fermion-spin framework, as a re- 
sult of self-consistent motion of the dressed holons and 
spinous, the incommensurate antiferromagnetism is de- 
veloped away from the half-filling, where the incommen- 
surate peaks are located at [(1 ± S)tt, tt] and [tt, (1 ± S)7r] 
with the incommensurability parameter S{x) defined as 
the deviation of the peak position from the AF wave 
vector [tt, tt]. This incommensurability parameter d{x) 
increases progressively with the doping concentration at 
lower dopings^^'^^, which leads to that the position of the 
lower-energy peak in k{uj) is doping dependent. Using 
a typical value of S in the underdoped regime obtained 
in the previous calculations^^'^^, the distance of the in- 
commensurate peaks is estimated as 20A. On the other 
hand, using a typical velocity of sound in La2-a;Srj;Cu04 
of about 5km/s and a phonon energy of 10 mev, one^^ 



4 



obtains a wavelength for the phonons of 15A. This value 
is in qualitative agreement with the distance of the in- 
commensurate peaks, and therefore the time scale of the 
dynamic incommensurate correlation is comparable to 
that of the lattice vibrations^^. In this case, the dy- 
namic lattice modulations are induced, then the dynamic 
spinon modulations dominate the heat transport of the 
underdoped cuprates, in other words, the incommen- 
surate spin fluctuation has been reflected in the ther- 
mal conductivity. On the other hand, these incommen- 
surate peaks are very sharp at low temperatures and 
energies^^'^^, however, they broaden and weaken in am- 
plitude as the energy increase for low energies w < xJ , 
and almost vanishes for high energies u > xJ, which 
leads to the shoulder appears in k{T) in the temperature 
T xJ. Finally, we emphasize that the present theory 
can describe the heat transport of the underdoped sin- 
gle layer cuprates, where only incommensurate spin fluc- 
tuation is observed^^. However, both incommensurate 
spin fluctuation and commensurate [tt, tt] resonance have 
been observed in the underdoped bilayer cuprates in the 
normal state^^ due to the bilayer splitting in the band 
structure^**. This resonance may lead some additional 
features in the thermal conductivity of the underdoped 
bilayer cuprates^^, and these and other related issues are 
under investigation now. 

In summary, we have studied the heat transport of the 
underdoped single layer cuprates within the t-J model. 
Our results show that the frequency dependence of the 
thermal conductivity spectrum k{uj) at low temperatures 
consists of two bands. The high energy band shows a 
weak peak, while the low energy peak is located at a 
finite energy. This high energy broad band is severely 
suppressed with increasing temperatures, and vanishes 
at higher temperature. Moreover, the temperature de- 
pendence of the thermal conductivity k{T) increases 
monotonously with increasing temperatures. Our results 
also show that although both dressed holons and spinous 
are responsible; for the thermal conductivity, the con- 
tribution from the dressed spinous dominates the heat 
transport of the underdoped cuprates. On the other 
hand, we emphasize that although the simplest t-J model 
can not be regarded as a complete model for the quantita- 
tive comparison with the doped cuprates, our present re- 
sults of the thermal conductivity are in qualitative agree- 
ment with the major experimental observations of the 
underdoped single layer cuprates^. 
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